O 



(1+1)-Dirac bound states in one-dimensional heterostructures; position-dependent 

Fermi velocity and mass 

Omar MustafcE 

Department of Physics, Eastern Mediterranean University, G. Magusa, north Cyprus, Mersin 10 - Turkey, 

Tel.: +90 392 6301378; fax: +90 3692 365 1604. 

Abstract: We extend Panella and Roy's 13] work for massless Dirac particles with position- 
dependent (PD) velocity. We consider Dirac particles where the mass and velocity are both position- 
dependent. Bound states in the continuum (BlC)-like and discrete bound state solutions are re- 
ported. It is observed that BIC-like solutions are not only feasible for the ultra-relativistic (massless) 
Dirac particles but also for Dirac particles with PD-mass and PD- velocity that satisfy the condition 
m {x) v% (x) = A, where > is constant. A Dirac Poschl- Teller and a harmonic oscillator models 
are also reported. 

PACS No.: 03.63. Fd, 03.65. Pm, 03.65.Ge 



, Keywords: (l-l-l)-Dirac equation, position-dependent mass and velocity, bound-states in the 

■ continuum, one-dimensional heterostructure 

o 

O: 
o 

I. INTRODUCTION 



] In heterostructure physics, it was believed that electrons are effectively described by the position-dependent mass 
O f (PDM) Schrodinger Hamiltonian (i.e., von Roos Hamiltonian, e.g., )■ Using the Pauli spin matrices in the 

Schrodinger Hamiltonian, the Dirac Hamiltonian was ignored. However, this perspective has drastically changed 
since the discovery of graphene [1, Q . Many studies on the applicability of Dirac Hamiltonian in condense matter 
^ ■ were carried out (cf., e.g., 4-1^ and related references cited therein) . It is found that the effective low-energy model 
, for the quasi-particles is ultrarelativistic (i.e., massless) and is described by the Hamiltonian 



H = VF(Tp. (1) 

■ Which is in fact the Dirac Hamiltonian for massless particles with an effective Fermi velocity vp (where vp = c/300, 

. c is the speed of hght, cr is a vector using Pauli matrices and p = — iV, with h—1). 

T-H ■ However, the information on the material properties may be encoded in the Fermi velocity of the Dirac particles 

I [1, 0|- In this case, the Dirac Hamiltonian (1) takes the form 



H = y^vp [x)a-j,'p-j,^VF {x). (2) 



■ Hereby, one should notice that the replacement of the constant velocity, wf, by the position-dependent one, vp (x)^ 
^ , would render Hamiltonian (1) non-Hermitian. Whereas the form of Hamiltonian (2) preserves Hermiticity and recovers 

■ the constant vp setting. Panella and Roy for example, have used Hamiltonian (2) to study bound states in the 

■ continuum (BIG) and discrete energy states for massless Dirac particle. Throughout this paper, we shall refer to their 
study as Panella-Roy's model (namely, their model with rn (x) = and vp (x) = vq cosh^ ax). They have found that 
with proper PD-Fermi velocity profile it is possible to create BIC-like and discrete bound-state solutions. 

In this paper, we extend Panella and Roy's [Tsj work and consider the Dirac Hamiltonian where the Fermi velocity 
and the mass are both position dependent. That is, we shall work with the Hamiltonian 

H = ^y Vp {x)axPx\/vF (x) + l3m (x) vp {xf , (3) 

where ax and /3 are the usual Pauli matrices The organization of this paper is in order. 

We discuss Hamiltonian (3) and give our methodical proposal in section II. We provide illustrative examples, 
including ultra-relativistic Dirac quasi-particles (i.e., particles with m(x) = 0), in section HI. In the same section, 
we show that similar scenarios (as those in the Panella-Roy's model [l^] for BIC-like and for discrete bound-states 
solutions) are observed for a wider class of m (a;) and vp (x) (i.e., for m{x)vp (x) = A, where ^4 > is constant). 
For such mass and Fermi velocity settings, a shift in the energy levels is obtained. In section IV, we show that 
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Dirac particles may be trapped in an effective Poschl- Teller potential [ij, ll5| produced by both their PD-mass and 
PD-Fermi velocity. Moreover, we show, in section V, that the (1+1)-Dirac oscillator may just be a consequence of a 
linear PD-Fermi velocity and a singular PD-mass (i.e., vp (x) — vqx and m (x) = A/x). Our concluding remarks are 
given in section VI. 



II. (l-M)-DIRAC PARTICLES WITH POSITION-DEPENDENT VELOCITY AND MASS 



With the usual textbook Pauli matrices and Dirac spinors, the (l-l-l)-Dirac equation Hip (x) — E'ip{x), for H in 
(3), would decouple into 



-iF{x)d.,[Fix)^^ix)] = Ci (2^)^1(2^) ; CAx)=E-mix)Fix)\ 
-tF{x)d^[Fix)i^,ix)] = (2^^2(2^) ; C2i^)=E + mix)F{x)\ 

where F (x) = \/vf {x) and 



'ij;{x)=N 



i>2 i^) J ' 



(4) 
(5) 



(6) 



(with N as the normalization constant) are used. Now let us multiply (4) and (5), from the left, by F (x) and use the 
substitutions ipi 2 (x) — F (x) ipi 2 (x) to imply 



^2 (x) -« 



C2 (•^) 



dx^i (x) . 



Which when substituted in (4) yields 



d3i (^) + 



C2 i^) 



- 2 



F'jx) 
F{x) 



d.j:^pi (x) = 



Ci (x) C2 (x) 



■01 {x) , 



(7) 



(8) 



where primes denote derivatives with respect to x. To get rid of the first order derivative and bring (8) into the 
one-dimensional form of Schrodinger equation we use 



V"! {x) = (2 i^Y (Pliq) ; q = q {x) . 

This would suggest that 

q' (x)=C2 {x)'-'- F {x)-\ 
However, one also needs to avoid position-dependent energies and choose = 1/2 to imply 

'3fC2{^)V lC2(a^) U^'Fi^)\C'2{^) 



(9) 
(10) 



df^b, {q)+vp {xY 



4VC2(a^)/ ^C^ix) 2\vf{x)J C^{x) 



with 



q{x) 



dy, 



^^{q)= E'-m{xYvF{xY cj^^iq), (11) 



(12) 



where q {x) represents a point canonical transformation. It is obvious that for massless particles equation (11) collapses 
into its most simplistic form 

-9^0^ (q) = E'^, (q) , 

that looks very much like the one-dimensional Schrodinger equation for free particles. However, the form of vp (x) 
would determine the domain of q (x) in (12) and, therefore, has its say in the process. This is to be clarified in the 
illustrative examples below. 



Nevertheless, one may use the non-relativistic limit m{x)vp{x) » Ei 



binding 



bind \ 



where E] 



'bind 



E 



m {x) Vp (x) (analogous to the textbook non-relativistic limit for Dirac particles with rest mass energy moC >> Eund 
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[Tg}). One would, in this way, recover the constant non-zero mass and constant velocity settings as well as accommodate 
position-dependent mass m (x) at vp {x) = c. This non-relativistic limit would, in turn, yield 

11 1 1 

(13) 



C2 i^) E + m {x) Vp (x)'^ Ebind + 2m (x) vp (x)'^ 2m {x) vp {xY 
Consequently, one may recast (11) as 



{Q) + Veff {q)4>i {q) = E^4>, {q) 



(14) 



where 



Veff (q) 



m (x) Vp (xY 



16 m (.t) Vp {xY 



(m{x)vp{xY'j fv'p{x)\ (m{x)vp{xf'j 



m (x) 



vp (x) 



m {x) 



+ m [xY V p {x^^ . (15) 



Obviously, this approximation may only be used for non-zero constant mass and not for massless Dirac particles. To 
illustrate our methodical proposal above we discuss the following illustrative examples. 



III. BIC-LIKE AND DISCRETE BOUND-STATES SOLUTIONS: PARALLEL AND COMPLEMENTARY 

TO PANELLA-ROY'S MODEL 

One consider the class of PD-mass and PD-Fcrmi velocity satisfying the condition m [x) vp {xY — ^, where A> 
is a constant. Under such assumptions, equation (11) would read 

- 9^01 (q) = {q) , (16) 

where A = E'^ — A^. Yet, one may notice that m (x) = is just a special case of the current more general proposal 
than that used in Panella- Roy's model ^3|]. Although this equation looks like Schrodinger equation for free particle, 
the domain of q (x) in (12) would determine the boundary conditions on the related state functions. This is to be 
clarified in the following two examples. The first of which is discussed here as a complementary model that reports 
on the consequences of using m (x) = mo/ cosh"* ax, and vp (x) = vq cosh^ ax settings on Panella-Roy's model. The 
second example considers vp {x) = vq (l + a^x^Y {x) — mo/(l -^a^x^Y shares similar scenario on the BIC-like 
and the discrete bound-states solutions as that reported in Panella-Roy's model 13]. 



A. Consequences of m{x)vF {xY = A on Panella-Roy's model: complementary 

Let the PD-mass and the PD-Fermi velocity take the following forms 

™ {x) — T — I and Vp {x) — Vq cosh^ ax (17) 

cosh ax 

respectively, with the constants mo, vq > and A = itIoVq. This would, in turn, imply 

r 1 1 

q{x) — / —dy= tanhax, (18) 

J VP (y) avo 

and suggest that q{x) € {—1/avQ, l/avo). Therefore, our particle under consideration is not free but rather quasi-free 
(i.e., trapped in a force field produced by its own PD-mass and PD-Fermi velocity in (17)) and is confined to move 
between —l/avo and +l/avQ. 

Whilst the solution of (16) is straightforward and takes the form 

0i(<z) =sinA<z, (19) 

it is rather unphysical (i.e., it does not satisfy the boundary conditions imposed by the range of q (x)). Nevertheless, 
one may use this unphysical solution to obtain the related un-normalized wave function components through (9) and 
(7) as 

; C2 = ^ + niovl (20) 



V-i (x) 



^1 (x) 

y/vYYxj 



-sech (ax) sin 
Wo 



■ tanh ax 



avo 
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and 



V'2 (x) ~ — — sech (ax) COS 

(a;) V vo 

Under such settings, the probabiUty density p (x) is given by 



■ tanh ( 



Ci^ E - moWg. 



2 2 C 

p{x) = \tPi{x)\ +\i^2ix)\ = —sech^ (ax) + 2moVQsech^ {ax) cos^ 

^0 



and the normalization constant N is obtained through 

t2 f Cl , ^UloV^ 



avQ 



■ tanh ax 



(21) 



(22) 



/ p (x) dx = 2N^(-^ 
J \avo 

— oo 

Moreover, the probabihty current density 



A 



= 1 



N = 



avQ\ 



\ 2 (aCi + 2amoV^ 



(23) 



jx = \/vF {x)'>j*i (x) y/vp {x)lp2 (x) + C.C. 



{i — i) sin ( tanh ax ] cos ( tanh ax 

\ avo / \ avQ 



= 0, 



(24) 



which is expected to vanish for bound states. As a result, the Dirac spinor in (6) and the related two components 
ipi [x) and 'tp2 (x) represent BIC-like solution. 

However, to make the unphysical solution in (19) satisfy the related boundary conditions (q) = at g (x) — 
±l/avo (hence becomes a physically admissible solution) one may shift q (x) — > q (x) + 1/avQ and recast the solution 
as 



01 (q) = sin 

This would immediately vanish at q = — 1/awo, and yield 



avoj 



An — 



mravQ 



En=± 



/mravQ \ ' 



; n = 1,2,3, ••• 



2 ' " V V 2 / 
for q = +l/avo. In this case, one obtains the un-normalized components as 



— — (tanh ax + 1) 
avo 



V'l (x) = \ — sech (ax) sin 



(25) 



(26) 



(27) 



-02 (x) = — sech (ax) cos 
V ^'o 



An 

avQ 



(tanh ax + 1) 



and the Dirac spinor would consequently read 

V'n (x) = A'n sech (ax) 



3in (tanh ax + 1)] 



^ cos (tanhax + 1)] 



(28) 



(29) 



Where A'n is given in (23). Yet, it should be noticed here that for mo = one may recover the final results of 
Panella- Roy's model fl3l| to obtain En = ±mravo/2, Nn = \J awo/2i?n, and 



V'n (x) 



a ( sin [if^ (tanhax + 1)] 

— sech (ax) 

2 cos [if (tanhax + 1)] 



(30) 
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As such, the what may look hke as a BIC-solution (documented in (19) and (24)) may turn out to be a bound 
state solution with discrete energy levels, if the proper physical boundary conditions are invested in the process 
(documented in (25) and (26)). Moreover, one observes a shift-up of order ttiIvq in the total energy squared, E^, 
and some scaling factors in the components of the Dirac spinor (i.e., yC^/^ "01 (^) S'^id \/Ci/wo for ip2 (x)), as 
discrepancies between our current model and Panella-Roy's model 13]. Obviously, should our m (x) = uio (i.e., the 
rest mass) and vp (x) — c (i.e., speed of light), then our q (x) — x/c and equation (11) would collapse into the regular 
textbook Dirac equation for a free particle where the total energy reads E — ±TOoC^. 



B. Parallel to Panella-Roy's model: vp (x) — vq [l + q^x^) 
We now consider that the PD-mass as 



m (x) 



(1 + a'^x'^y 



and the PD- Fermi velocity as 



Vp {x) = Uo (l + a^x^) q {x) 



avo 



■ arctanax. 



(31) 



It is easy to observe similar scenario as that associated with (j>i (q) of (19), where in the current case the particle 
described in (16) is now confined to move within — 7r/2az;o ^ 9 (x) < t: /2avQ. The unphysical solution then reads 



^1 (9) — ^9 — { arctan ax 



(32) 



Which would, in turn, imply that 



V-i (x) 



^i(^) _ fc 



\Jvf {x) V ^'0 Vl + oP-x 



: sm 



■ arctan ax 



avo 



and 



Therefore, 

j Pix)dx 

— CO 

and 



1P2 {x) 



^2 {x) 

\/vF [x) 



Ci 1 



arctan ax 



vo (1 + a'^x'^) 



(^i + 2moVQ sin^ ( arctan ax 



avo 



dx = l=^ N = 



(33) 



(34) 



avQ 



(35) 



{i — i) sin arctan ax cos arctan ax 

\ avo / \ avo , 



0, 



(36) 



indicating the existence of bound states. As such, the Dirac spinor in (6) and the related components ip-^ (x) and 
tp2 (x) represent a BIC-like solution. 

However, the physically admissible solution would be achieved through a shift in q (x) — > q (x) + TT/2avQ to read 



^ (9) = sin 



and yields 



Xn = navo En = ±y (navo) + mlvQ ; n = 1, 2, 3 



(37) 



(38) 
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Therefore, 



and 



-02 (x) 



1 



: sm 



( arctan ( 



Consequently, the Dirac spinor would read 



n ^arctan ax + — ^ 



1 



V 



n = 1,2,3, ■ 



^ -^/l^sin [n (arctan ax + ^)] 
cos \_n (arctan ax + ^)] 

y^sin [n (arctan ax + §)] 



Moreover, for the case when m (x) = one may obtain = ±a/ {navo) and 



1 / sin [n (arctan ax + -I )] 
TT Vl + Y -i cos [n (arctan ax + |)] 



(39) 



(40) 



(41) 



(42) 



Again one observes similar effects of the m (x) vp {x) = A setting on the total energy and on the components of 
the Dirac spinor as those mentioned in the above example. 



IV. (1+1)-DIRAC POSCHL-TELLER HOLES FOR vf (x) = vo AND A PDM m{x) ^ 
Let us now consider the case where vp (x) = uq => q (x) = x/uo and 



m (x) 



sm ax 

Under such settings, the effective potential in (15) would read 



Veff {q) 



a 



(43) 



(44) 



Which 
(14) as 



16 sin^ (avoq) 

is obviously a shifted Poschl- Teller type periodical potential (cf., e.g., [13, [l^)- In this case, one may rewrite 



Vo 



sis-l) 



sin^ (avoq) 



K (q) =[E'- ^) 01 (q) 



(45) 



where Vo = 2a'^VQ and s(s — 1) = m'^/a^ — 5/16. Such periodical potential setting imposes infinite impenetrable 
barriers manifested by the singularities between the holes (i.e., at g = 0, 7r/woa, 27r/uoa, • • • or equivalently at x = 
0, 7r/a, 27r/a, •••.). Here we pick up the hole within < = vaq < 7r/a to obtain 



01 (g) = sin"* {avnq) 2F1 ( -n, s + n, s + ^; sin^ {av^q) 



(46) 



and 



^„ = ±^^l + 16(s + 2n)2 ; n = 0,l,2,. 



(47) 
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Where 



s^^ + ^mla^- (48) 



Then, one would, in a straightforward manner, cast 



V'l {x) = Nn\ — + sin-' (ax) 2F1 ( -n, s + n,s+ ^; sin^ (ax)] (49) 

V vq smax \ 2 J 

and find V'2 (2^)? using (7), to construct the Dirac spinor of (5). Obviously, BIC-like bound states are not feasible here 
and only discrete bound state solutions are obtained. 



V. (1+1)-DIRAC EFFECTIVE HARMONIC OSCILLATOR TOY: A BY-PRODUCT OF vf (x) = vox AND 

A m{x) = A/x 

Consider a singular position-dependent mass along with a linear position-dependent Fermi velocity of the forms 



m{x) — A/x; vp (x) = vqx ; x € (O.oo) = 
In this case, equation (14) along with (15) would yield 

1 



' (x) = lnx-^/^°; q e (-00, 00) . 



d'^4>i (q) + -yq'<Pi (g) = ( i?^ + ^ ) 01 {q) 



(50) 



(51) 



with w = 2Avq. Obviously, 4>i {q = zLoo) = represent the boundary conditions for the current Dirac harmonic 
oscillator at hand. In a straightforward manner, one would use the traditional textbook procedure and find that 



16 



= uj I n + 



1 



En=±\lAvl{2n + l)-^. 



and 



where _ff„ (^y^ Av^q^ are the Hermite polynomials. Then we may obtain 



" VqX 



(52) 



(53) 



(54) 



and find "02 {^) using (7) to construct the Dirac spinor of (5). Only discrete bound state solutions are observed here. 



VI. CONCLUDING REMARKS 



In this work, we have considered the (1+1)-Dirac particles where the mass and the Fermi velocity are both position- 
dependent. An alternative methodical proposal is proposed in such a way that the Panella- Roy's model [l^ becomes 
a special case. The set of m {x) and vp {x) that satisfies m{x)vF {x) = A is a wider set than that used by Panella 
and Roy who have used massless Dirac particles. Moreover, analogous to the well known textbook non-relativistic 
limit for Dirac particles (i.e., rest mass energy moC^ >> E^ind, where Ebind = E — irioC^), we have used the limit 
where m (x) vp {x)^ » Eund for non-zero PD-masses. To the best of our knowledge such methodical proposal has 
not been reported elsewhere. 

For Dirac particles with m{x) and vf{x) satisfying m{x)vp{x) = A, we have reported feasible BIC-like and 
discrete bound-states solutions (documented in section III). They are in an almost exact accord with the scenario 
reported in the Panella-Roy's model. However, we have also observed a shift-up of order m^WQ = A^ in the total 
energy squared, E"^, and some scaling factors in the components of the Dirac spinor (i.e., y/Cih^ for ijji (x) and 

Ci/^^o for 0^2 {^))- Moreover, the results of our methodical proposal collapse into those of Panella and Roy in [l3| 
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for TTio = 0. Yet, should one use m (x) = mo (i.e., the rest mass) and vp (x) = c (i.e., speed of hght), then q (x) = x/c 
and equation (11) would collapse into the regular textbook Dirac equation for free particle, where the total energy is 

E — ±TOoC^. 

Finally, for the case where m (x) vp {x)"^ A, we have shown that Dirac particles may be trapped in an effective 
force fields produced by both their PD-mass and PD-Fermi velocity. This is documented in the effective Poschl- Teller 
and the effective harmonic oscillator models discussed in sections IV and V, respectively. No BIC-like bound state 
solutions are observed for these models. 
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